Quantum localization in rough billiards 
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We study the level spacing statistics p{s) and eigenfunction properties in a billiard with a rough 
boundary. Quantum effects lead to localization of classical diffusion in the angular momentum space 
and the Shnirelman peak in p{s) at small s. The ergodic regime with Wigner-Dyson statistics is 
identified as a function of roughness. Applications for the Q-spoiling in optical resonators are also 
discussed. 
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PACS numbers: 05.45. +b, 72.15.Rn, 42.55.Mv 

In 1984, Bohigas, Giannoni, and Schmit J^] demon- 
strated that random matrix theory (RMT) [^| describes 
the level spacing statistics of classically chaotic billiards. 
After that such types of billiards have been studied in 
great detail by different groups However, all billiards 
under investigation were characterized only by one typi- 
cal time scale tc, namely, the time between collisions with 
the boundary. Another type of chaotic systems with dif- 
fusive behavior, like the kicked rotator Q|, has an addi- 
tional much longer time scale to ^ tc which is required to 
cover the accessible phase space. In this situation quan- 
tum interference effects may lead to exponential dynam- 
ical localization of the eigenstates and disappearance of 
level repulsion. 

Recently, it has been shown ||] that the diffusive 
regime also appears in a nearly circular Bunimovich sta- 
dium billiard. The authors of [|| gave an estimate for 
the localization length in the angular momentum space 
and found the energy border E^rg above which the eigen- 
states become ergodic on the energy surface. Their nu- 
merical results demonstrate the change of level statistics 
from Wigner-Dyson to Poisson when the energy becomes 
smaller than E^rg- However, this example being very 
interesting for mathematical studies is not physically re- 
alistic. 

At the same time a great progress has been reached in 
optics of microcavities like micrometer-size droplets ^] 
and microdisk lasers (t). The industrial request to pro- 
duce directed light pushed the researches to investigate 
ray dynamics in weakly deformed circular billiards and 
droplets. It was shown that above some critical deforma- 
tion the ray dynamics becomes chaotic. As a consequence 
the directionality of light from droplets and Q-factors of 
such resonators are strongly affected However, due to 
smoothness of the deformation the diffusive regime was 
hardly accessible in such systems. 

In this Letter, we investigate another type of generic 
boundary deformation which may have important phys- 
ical applications. Namely, we consider billiards with a 
rough boundary. In this approach, the boundary is a 
random surface with some finite correlation length. The 
physical realizations of such situation can be quite dif- 



ferent. As examples, we can mention surface waves in 
the droplets which are practically static for the light [D, 
nonideal surfaces in microdisk lasers 0], and capillary 
waves on a surface of small metallic clusters On a 
first glance it seems that such a rough boundary in a 
circular billiard will destroy the conservation of angular 
momentum leading to ergodic eigenstates and RMT level 
statistics. In spite of this we show that there is a region of 
roughness in which the classical dynamics is chaotic but 
the eigenstates are localized and the level spacing statis- 
tics p{s) has the sharp Shnirelman peak at small spacings 
s ||l^,|ll| . We also demonstrate the close relation between 
this model and the kicked rotator. 

As a model of a rough billiard we chose a circle with a 
deformed elastic boundary given by R{9) = Rq + AR{9) 
with AR{9)/Rq = Re J2m=2 e'^"^^ and where j„i are 
random complex coefficients. The roughness of the sur- 
face is given by k{6) = {dR/d9)/Ro. In the following, 
we will consider the case of weak roughness k <C 1. Fur- 
thermore, one can model different type of surfaces by 
chosing the appropriate dependence of the amplitudes 
on m. However, our results show that in the domain of 
strong chaos the classical diffusion and quantum localiza- 
tion in orbital momentum space are determined by the 
angle average = {K^{9))e- Due to that we choose a 
typical dependence 7^ ^ 1/m such that all harmonics 
give the same contribution in the roughness. 

First, we consider the classical ray dynamics which for 
K <C 1 can be described by the following rough map: 



1 = 1 + 2^1 



-I- TT - 2m:CS\l\{l/lmax) 



(1) 



Here the first equation gives the change of the angular 
momentum (and velocity vector) due to the collision with 
the boundary and the second one the change of angle be- 
tween the collisions. This map describes the dynamics 
in the vicinity of the resonant value 1^ defined by the 
condition 9 = 9-\-2'Kr with integer r. Our numerical sim- 
ulations of the exact ray dynamics show that the rough 
map (^) indeed gives an excellent description (see Fig. 
|l|). A similar map for a stadium billiard was given in 
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However, in contrast to where k{9) has a discontinu- 
ity, the global chaos sets in only if the roughness is above 
some critical value ^/MAR/R ^ k > Kc- Below Kc the 
KAM theory is valid and the phase space is divided by 
invariant curves. The chaos border can be estimated on 
the basis of Chirikov criteria of overlapping resonances 
[ p2[ which gives Hc ^ 4:M~^/'^ (the numerical coefficient 
is extracted from the data for M = 20). This border 
drops strongly with M and therefore we will concentrate 
on the regime of strong chaos without visible islands of 
stability (case of Fig. |l|). In this regime the spreading 
in angular momentum space goes in a diffusive way with 
a diffusion constant that is easily estimated from (Q) as 
D = [Alf/At = 4 {ll^^^ - ll) k2 where the time t is 
measured in number of collisions. 
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FIG. 1. Comparison of the exact dynamics (points) and 
the rough map (|^) (full curve k,{6)). The points correspond 
to a; = (F- l)/{2Vl'Lax - P) and M = 20, k = 0.011. 

Now, we turn to the investigation of the quantum prob- 
lem. For this we expand the wave function ipir, 9) with 
energy E — h^k^/2m in terms of the Hankel functions 
which form a complete set: 

^(^^ ^) = E Hl+\kr) e^'» -I- bi Hl-\kr) e''«) . (2) 

The regularity of -0 at r = requires a; — bi so that 
the ai are the amplitudes in angular momentum space. 
The boundary condition tlj{R{9),9) — results in a sec- 
ond equation hi = ^H' i^) o-i' where Sw {E) has the 
meaning of a scattering matrix for waves reflected at the 
rough boundary. The energy eigenvalues are determined 
by det[l — S{E)] = 0. A convenient expression for the 
S-matrix can be obtained by the quasiclassical approxi- 
mation H^^\kr) w 2 [27rfc;(r)r]-i/2 exp[±?(^/(r) - 7r/4)] 
where ki{r) — k{\ — rfjr^Y^'^i /^/(^) = ll^ dfkiir) and 
Ti « \l\lk is the classical turning point. From this rep- 
resentation and the boundary condition for i]} we obtain 
(for more details see [|3|): 

5'^^, ^ giMi(-Ro)+Vi'(-Ro)-»'r/2 ^ ^|gi2/c,^(_Ro)A_R(e)|^/ ^ 

(3) 



This is a local unitary expression for S near l^. which 
is valid for A? = |Z — <C Imax = kRo and M < Imax- 
A stationary phase approximation for the ^-integral gives 
the classical change of I [see (||)] and determines the struc- 
ture of S'-matrix. In fact, this matrix is very similar to the 
evolution operator of the kicked rotator Q corresponding 
to Ai? (X cosO. According to this analogy the localiza- 
tion length t is determined by the classical diffusion rate 
I — (3D/2 where (3 is the symmetry index for orthogo- 
nal (/3 = 1) ijl^ or local unitary symmetry {(3 — 2) [ p5[ . 
Since generally AR{9) ^ AR{-9), we have = 2 so 
that the localization length is directly determined by the 
roughness 

e = D = A{ll^^-ll)k^ (4) 

This result can also be derived on a more rigorous ground 
based on the supersymmetry approach for a model with 
random phases /i/ [ p^ . The expression (^ is only valid 
for D > M while 1 < D < AI corresponds to a 
more complicated regime with £ ^ M (see below). For 
1 < Z) <C Imax the eigenfunctions are exponentially lo- 
calized in the orbital space |a;| cx exp(— |Z — lo\/i) while 
the classical dynamics is ergodic on the whole energy sur- 
face. The quantum dynamics becomes ergodic only for 
^ > imax- According to the Weyl formula the level num- 
ber at energy E is N 'mRlE/2h? = Imaxl^- Therefore 
the states are ergodic for 



This border is much higher than the perturbative border 
N < Np « 1/(16k^) where the diffusion mixes less than 
one state (Z? w 1). These different regimes are presented 
in Fig. |. 
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FIG. 2. Global diagram showing the eigenstate properties 
for different values of level number A*' and roughness k for 
D 3> M. The two full lines give the ergodic (A^e) and per- 
turbative (Np) borders, the dashed line is the classical chaos 
border Kc « 0.002 for M = 20. The parameters of Figs. ^ 
are shown by the two points with arrows. 
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To check the above theoretical predictions, we have 
solved numerically the boundary condition for ■(/'(r, 9) 
(^. In this way, we have obtained both the energy 
eigenvalues and the amplitudes a; (with normalization 
J2i kip = !)• In the ergodic regime N > N^, we find 
that the level spacing statistics p{s) is in a good agree- 
ment with RMT (see Fig. ^) . On the contrary, in the lo- 
calized regime N < TVe, approximately each second level 
is quasidegenerate leading to the Shnirelman peak 
at small spacings (Fig. ||). This peak represents ap- 
proximately a fraction a « 0.33 of all spacings. The 
other spacings are described by a rescaled Poisson distri- 
bution p(s) = (1 — a)^ exp[— (1 — a) s]. The appearance of 
the Shnirelman peak is in agreement with the prediction 
made in . Its physical origin is the time reversal sym- 
metry {S-i^-i' = S*i, or a-i = a*) due to which two states 
localized around Ig and —Iq form a quasidegenerate pair 
of symmetric and antisymmetric states . The original 
Shnirelman theorem was formulated for quasi-integrable 
billiards Our case corresponds to the chaotic do- 

main, however, due to localization the peak still exists. 
A similar situation for the kicked rotator was studied in 
PH . The fraction of non degenerate levels (1 — 2a ~ 0.34) 
is due to states localized near Iq w 0. The measure of such 
states is approximately M/lmax ~ 1 — 2a where the nu- 
merical coefficient was extracted from our data (see Fig. 

This peak was not found in ||^ because the stadium 
billiard has additional symmetries and only the states of 
one parity were considered. 




FIG. 3. Level spacing distribution p{s) for Imax ~ 100, 
M = 20 and average roughness: a) k ~ 0.16 (ergodic, 
i ^ 1000) and b) k « 0.012 (localized, ^ ^ 6). The to- 
tal statistics is 2000 levels (2150 < iV < 2350) from 10 dif- 
ferent realizations of the rough boundary. Also shown are 
the Wigner-Dyson and rescaled Poisson distributions with 
a = 0.33. 

In Fig. ^, we show three typical eigenfunctions in an- 
gular momentum representation. In the ergodic case, 
the probability is homogeneously distributed nearly in 
the whole interval {—Imax, Imax)- We mention that in 
our computations we have included about 10 evanescent 
modes (with \l\ > kRo). The two other wave functions 



are in the localized regime of Fig. g with £ ^ Imax- One 
of them has a double peak structure due to tunneling 
between time reversed angular momentum states. This 
state corresponds to a quasidegenerate Snirelman state. 
We also observed many other similar states with much 
smaller level splitting (s < 10^^). The other localized 
state has the maximum « and corresponds to a non- 
degenerate level. Both states clearly show exponential 
localization indicated by the dashed lines. The numeri- 
cal values of the localization length are a bit higher than 
the estimate (01) with Ir ^ Iq Imax- 




FIG. 4. The amplitudes \ ai \ for three typical eigenfunctions 
corresponding to the parameters of Fig. H (with level number 
N ~ 2250): a) ergodic case; b) localized quasidegenerate state 
with s = 0.1 and fitted localization length £ ~ 9 (dashed line); 
c) nondegenerate localized state with fitted £ ~ 7.5 (dashed 
line). The full line shows the theoretical length [see Eq. (^] 
1^5. The rough boundary realization is the same as in Fig. 
|l (b,c) or appropriately rescaled (a). The curves b and c are 
shifted by factors 10~^ or 10~^, respectively. 
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FIG. 5. Rescaled localization length £/M as a function 
of rescaled diffusion rate D/M for Imax = 100: crosses 
are numerical data for billiards with 0.008 < k < 0.027 
and 20 < M < 60; squares are data for the effective 
kicked rotator model (see text) with 0.006 < k < 0.043, 
5 < Af < 60 and a matrix size of 600; the full curve is the 
fit ^ = 0.5Af/ln[l + 0.5/(0.05 + D/M)] with the asymptotic 
limit £ = D shown by the dashed line. 



3 



To understand this discrepancy, we have calculated i 
for a wide region of M and k. Many cases were sys- 
tematically different from (^) . To increase the parameter 
range we also studied the effective kicked rotator model 
with random phases iii in the evolution operator (and 
^max — 100, Ir = 0). The rcsults of these studies demon- 
strate the scaling behavior IjM = f{D/M) (see Fig. H). 
We attribute the remaining difference between the two 
models to a smaller sample size for the billiard (100 vs. 
600) and finite Ir ^ Imax values there. For D ^ M the 
scaling reproduces the estimate (||) while for D < M the 
length £ remains close to M. In this case, the evolution 
operator becomes a band random matrix of width M 
with strong diagonal fluctuations. They are much larger 
than the off-diagonal matrix elements so that £ ^ M for 
1< D< Af 0. 

Let us now discuss how the above properties of eigen- 
states will affect the characteristics of resonators with 
rough boundaries. In optical resonators the rays with 
large reflection angle escape from the system (see for ex- 
ample (H). Due to that there is an effective absorption 
in the momentum space for I < Ic where the critical mo- 
mentum is determined by the index of reflection so that 
typically Ic/lmax ~ 1/2. This absorption affects the Q- 
value of the resonator which is approximately equal to 
the number of collisions until the escape. In the ergodic 
regime this number is determined by diffusive spreading 
in the momentum space so that Q ~ Ic/D ~ being 
proportional to the inverse Thouless energy. On the other 
hand in the localized regime the probability to reach Ic is 
exponentially suppressed. The subsequent estimate gives 
InQ ~ lc/£ ^ i/{lmax fi"^) > 1 and determines the rough- 
ness border kq in the Q-spoiling which corresponds to 
our ergodic border (^. For given Imax this border is 

i^Q ~ 1/ (2 \^lmax)- 

Above we discussed the effects of the rough boundary 
in a circular billiard. However, we can argue that sim- 
ilar effects should be observed in more general types of 
convex smooth billiards. Indeed, in such systems, a large 
fraction of the phase space is integrable and character- 
ized by two quantum numbers. As in the circular case 
the rough boundary will lead to a diffusive behavior in 
one of the quantum numbers parallel to the energy sur- 
face and again the quantum interference effects can give 
its localization. We note that our results are valid for 
weakly rough billiards k <C 1 while the case of strong 
roughness k ^ 1 deserves separate studies. 

A generalization to the three dimensional (3d) case 
represents an interesting direction for further research. 
We expect that for a sphere with rough boundary the 
above 2d analysis is very relevant. Indeed, in a perfect 
sphere a trajectory is confined to a plane and the preces- 
sion frequency of this plane is zero. Due to roughness this 
plane will slowly precess with a frequency proportional to 
K. This adiabatic process will weakly affect the dynamics 
and quantum localization inside the plane section. The 



localization should give rise to the Shnirelman peak in 
3d. However, the situation should be quite different in 
more general 3d smooth integrable billiards (e.g. ellip- 
soid). There the precession frequency is rather high and 
the diffusion becomes really two dimensional on the en- 
ergy surface. As for locahzation in 2d disordered systems, 
one might expect an exponential fast increase of £ with 
the roughness {\n£ ~ /^^^ k^). 

In conclusion, we studied the billiards with strongly 
different diffusive and instability time scales. In this sit- 
uation quantum effects can produce localized eigenstates 
breaking classical ergodicity on the energy surface. We 
wonder if the periodic orbit approach can lead to a deeper 
understanding of the physical situation. 
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